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O 1 INTRODUCTION 



ABSTRACT 

Hoyle & Folwler showed that there could be Radiation Pressure Supported Stars 
(RPSS) even in Newtonian gravity. Much later, Mitra found that one could also con- 
ceive of their General Relativistic (GR) version, "Relativistic Radiation Pressure Sup- 
ported Stars" (RRPSSs). While RPSSs have z < 1, RRPSSs have z > 1, where z 
is the surface gravitational redshift. Here we elaborate on the formation of RRPSSs 
during continued gravitational collapse by recalling that a contracting massive star 
must start trapping radiation as it would enter its photon sphere. It is found that, 
irrespective of the details of the contraction process, the trapped radiation fltix should 
attain the corresponding Eddington value at sufBciently large z ^ 1. This means that 
continued GR collapse may generate an intermediate RRPSS with z ^ 1 before a 
true BH state with z = oo is formed asymptotically. An exciting consequence of this 
is that the stellar mass black hole candidates, at present epoch, should be hot balls 
of quark gluon plasma, as has been discussed by Royzen in a recent article entitled 
"QCD against black holes?\ 
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1 For any self-gravitating object, the definition of "compact- 
. 5^ ' ness" may be given in terms of the surface gravitational 
^ . redshift (Weinberg 1972): 



z = {l-Rs/R) 



-1/2 



(1) 



where Ra = 2GM/(? is the Schwarzschild radius of the ob- 
ject having a gravitational mass M and radius R. Here, G 
is the Newtonian gravitational constant and c is the speed 
of light. For the Sun, one has z ~2 x 10~^ while for a typ- 
ical neutron star z ~ 0.15. Another important parameter 
for a self-gravitating object is a; = Pr/Pg, where pr is the 
pure radiation pressure and Pg is the kinetic pressure (Mi- 
tra 2009a,b). While the central region of Sun has x ~ 0.006, 
obviously, a cold object at temperature T = Q has x = Q. 
It is known that strictly static and cold stars have an upper 
mass limit both in Newtonian and Einstein gravity. In the 
Newtonian case, this is obtained by the marriage of Newto- 
nian gravity and Special Relativity, and is known as "Chan- 
drasekhar Mass", M^h- On the other hand, in Einstein grav- 
ity, such an upper mass limit is known as "Oppenheimer - 
Volkoff Mass" (Weinberg 1972). The precise values of such 
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limits depend on the equation of state (EOS) and other de- 
tails. However, for a pure cold, Helium dwarf, M^h ~ 1.4M©, 
where Mq is the solar mass. On the other hand, for a pure 
free neutron Fermi-Dirac fluid, Mov ~ O.SM©. But if such 
upper limits were the full story, we would not have had stars 
of masses as large as ~ IQQMq. Further, there are inter- 
stellar gas clouds of mass probably as large as ~ 10® Mq. 
The reason behind the existence of non-singular cosmic ob- 
jects with such higher masses is that they are not supported 
by cold quantum pressure alone. On the other hand, they 
are supported not only by pg but partly by pr too (Mitra 
2009a,b). 

Even when one would work at a purely Newtonian level 
{z <^ 1), the probable increase in the value of x would sup- 
port the higher self-gravity of a star M 3> IMq. And this 
was probably first realized by Hoyle & Fowler (1963) and 
Fowler (1966) who conceived of the Radiation Pressure Sup- 
ported Stars (RPSS) having a; ^ 1. It turns out that, 
given the self-imposed restriction, z <^1, one would require 
M > 72OOM0 in order to have x>l (Weinberg 1972). With 
the increase of x and attendant self-gravity, the star tends 
to be more compact, i.e., z would increase Mitra(2009a,b). 
For instance even though "supermassive stars" are Newto- 
nian (i.e., z <C 1), they possess z much larger than the solar 
value of ~ 2 X 10~®. Accordingly, it is possible to con- 
ceive of a Newtonian RPSS with z ~ 0.1, i.e., one which is 
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almost as compact as a neutron star. A related important 
concept here is the "Eddington Luminosity" . In order that a 
self-gravitating and self -luminous object can remain quasi- 
static, its luminosity at a given radius r = r must be less 
than a critical value (Hoyle & Fowler 1963, Weinberg 1972): 



Led{r) 



47rGM(r)c 



(2) 



where M(r) is the gravitational mass within a given surface, 
and K is the appropriate opacity. In case, the star will be 
composed of pure ionized hydrogen, one will have the lowest 
value of opacity called Thomson opacity: k = ar/nip ~ 0.4 
cm^ g~^ where, ar is the Thomson crosssection and nip is 
the proton rest mass. If one would have the intrinsic lumi- 
nosity, L > Led, the star would be disrupted by radiation 
pressure. The Newtonian supermassive stars, conceived by 
Hoyle & Fowler necessarily radiate at the above mentioned 
Newtonian Eddington rate. As one would tread into Einstein 
gravity, however, one should be able to conceive of situations 
with 2^1. The definition of Eddington luminosity, in such 
a case, gets modified (Mitra 1998): 

Led = (1 + ^) (3) 

Accordingly, in Einstein gravity, there should be Rclativis- 
tic RPSSs (RRPSS) which would locally radiate at the 
above mentioned enhanced Edddington rate. In the follow- 
ing, we discuss the basic physics behind the formation of 
such RRPSSs. 



2 MORE ON EDDINGTON LUMINOSITY 

It is known that the concept of Eddington luminosity is rele- 
vant not only for the structure of stars but for the accretion 
process around the stars too. In fact, the accretion luminos- 
ity of the star also must be limited by L^d for steady spher- 
ical accretion. To appreciate this, let us recall here the basic 
physics behind the concept of "Eddington Luminosity" by 
considering the fluid to be a fully ionized hydrogen plasma: 
The average attractive gravitational force on one atom is 
Fg = —GMnip/FI?. If this would be the ordy force acting 
on the plasma, then both the intrinsic luminosity and the 
accretion luminosity could be infinite. In reality, however, L 
is finite and has a maximal value, Led, because the ionized 
H-atom is also subject to the repulsive force due to accretion 
luminosity of the central object: 



Frad = or q/c 



(4) 



where q = L/AitB? is the radial energy /heat flux. Thus the 
effective, 

Fg^Fg+Frad = -^{{l-OL)mp\ (5) 

where a = L/Led and Led = ^nGMc/ Therefore, as if, the 
radially outward heat transport in a spherically symmetric 
isotropic isolated body reduces the Effective Gravity (EG) 
by a factor of (1 — a). This may be also seen as a reduction of 
the inertial mass (IM) by the same factor. A value of a > 1, 
i.e., L > Led would thus mean EG to bo negative, where 
either the star would be disrupted despite self-gravity or 
instead of accretion there could be radiation driven winds. 
The latter indeed happens in the atmosphere of very massive 



stars with very largo L. In the context of collapse, if the 
collapse generated luminosity would approach its Eddington 
value, it becomes clear then that the collapse process would 
tend to be stalled. Essentially, Led corresponds to a critical 
comoving outward heat flux of 



9ed = 



GM 



(1 + ^) 



(6) 



47r7?2 i^R2 

Irrespective of the specific mode of reduction of the EG due 
to heat flow, the very notion of an "Eddington Luminos- 
ity" , both in Newtonian gravity and in GR, implies that the 

attainment of L — Led would stop inflow/ collapse. In fact, 
in a very important study on GR collapse, Herrera & San- 
tos(2004) have indeed shown that outward heat flow reduces 
IM by a factor (1 — q). Therefore, in principle, the GR col- 
lapse process can certainly slow down and get stalled if the 
collapse generated luminosity would approach its maximal 
value (a — > 1)! Further numerical as well as analytical stud- 
ies of radiative GR collapse have confirmed the above men- 
tioned analytical result (Herrera & Santos 2004; Herrera, 
Frisco & Barrcto 2006; Herrera, Frisco, & Ospino, 2006; Her- 
rera, Frisco, Fucnmayor, & Troconis 2009). If such papers 
(Herrera & Santos 2004; Herrera, Frisco & Barreto 2006) 
would be used for weak gravity, one would clearly identify 
the Q occurring in them as none other than the a appear- 
ing in Eq.(5). While such studies considered enhancement of 
radiation flux due to matter -radiation interaction and are 
somewhat non-generic, we would consider here a generic ef- 
fect: the unabated enhancement of the gravitationally trapped 
radiation beyond the photon sphere and consequent attain- 
ment of Eddington luminosity. 



3 SELF-GRAVITATIONAL TRAPPING OF 
RADIATION 

General Relativity (GR) predicts that even the trajectories 
of quanta emitted by a star itself do bend away from the 
direction of normal to the direction of the tangent of the 
surface of the body because of the effect of the gravitational 
field of the star. However, as long as 2 < yfi— 1, the emitted 
quanta nevertheless manage to evade entrapment and move 
away to infinity. But if the body would be so compact as 
to lie within its "photon sphere", i.e., R < (3/2) Jis or ^; > 
\/3— 1, then only the radiation emitted within a cone defined 
by a semi-angle 9c (Harrison 2000): 



{l-Rs/Rf^\Rs/R) 



(7) 



will be able to escape. Radiation emitted in the rest of the 
hemisphere would eventually return within the compact ob- 
ject. 

In the presence of an external radiation, i.e., in the 
absence of a strict exterior vacuum, even an apparently 
static supermassive star is not strictly static. This is so bo- 
cause, the spacetime in such a case is described by radia- 
tive Vaidya solution (Vaidya 1951) rather than the exact 
vacuum Schwarzschild solution. In this sense, while the ex- 
terior spacetime of a cold White Dwarf or a Neutron Star 
having a fluid at zero tompetaturc is described by vacuum 
Schwarzschild metric, the exterior spacetime of any radia- 
tive object including the Sun or a supermassive star, in a 
strict sense is described by the Vaidya metric. Accordingly, 
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in a strict GR sense, stars are in quasistatic equilibrium and 
always evolving. For the Vaidya metric, the value of z could 
be arbitrary high, and that is the reason a collapsing and 
radiating object is supposed to attain the black hole (BH) 
stage having 2 — 00. As the collapse/contraction proceeds, 
both M and R decrease so that z increases and eventually, 
the collapsing body must approach the exact BH state with 
z = 00. The effect of radiation trapping during these in- 
termediate high-2 states has never been considered, though 
previously Kembhavi & Vishveshwara(1980) observed that: 

"If neutrinos are trapped, they will not be able to trans- 
port energy to the outside, and this can have serious conse- 
quences on the thermal evolution of the star. These consid- 
erations might become especially interesting in the case of a 
collapsing phase which leads to the formation of a compact, 
dense object." 

At high z, J? « i?s and from Eq.(7), one can see that, 
sin6'c 6c ~ (\/27/2)(l -|- z)~^. Hence the solid angle of 
escaping radiation is 

.2 277r , 



-(1 + ^) 



(8) 



The chance of escape of radiation therefore decreases as 
Qc/^n « (27/8)(l -I- z)-^. This means that if without trap- 
ping 10^" neutrinos/photons would escape a particular spot 
on the surface, with gravitational trapping, only 1 out the 
10^" quanta would escape for z = 10" . Consequently, as the 
collapse generates internal heat/radiation, the energy den- 
sity of trapped radiation pr and associated outward heat 
flux within the body would increase as 

(9) 



Qtra 



R-^il + zf 



The distantly observed luminosity /flux would be lesser by 
a factor of (1 -|- 2)^ than what is indicated by the foregoing 
Eq. Such a reduction would take into account the fact that, 
locally, trapped quanta are moving in almost closed orbits. 
But as far as local flux is concerned, to avoid double count- 
ing, one must not introduce any additional factor of (1 -f- 2) 
in Eq.(9). Using Eqs.(6) and(9), we see that, in this regime 
of 2 > 1, 

g*^^(l_f2) 

Qed RM ^ ' 

Initially, of course, a <C 1. But during the collapse, both 
R and M would decrease and Eq.(lO) would show that, as 
2 — ^ 00, a would increase dramatically, a 1, at a suffi- 
ciently high 2. At this stage the collapse would degenerate 
into a secular quasistatic contraction by the very definition 
of Led- As if, a leaking and contracting balloon stops con- 
traction as its self gravity fixes the leakage by forcing the 
molecules to move in (almost) closed circular orbits. Also, si- 
multaneously, the attendant heat and pressure become large 
enough to resist further contraction. In a very strict sense, 
however, the body would still be contracting on extremely 
long time-scales! This is so because as long as an horizon 
is not formed, i.e., 2 < 00, the body would radiate and M 
would continue to decrease. Consequently the metric would 
remain non-static and, in response, R too, would decrease. 
It is this infinitesimal decrease in the value of R and atten- 
dant much higher secular increase in the value of 2 and qtrap 
which would generate just enough energy (at the expense of 
Mc^) to maintain the GR Eddington luminosity seen by a 
distant observer: 



■•ed 



AwR^q^d _ AnGMc 

(1-1-2)2 ~ k{1 + z) 



(11) 



Since L°° = —c^dM/du, the observed time scale associated 
with this phase is 



KC(1 -I- 2) 



-c^dM/du 47rG 



(12) 



Obviously, u 00 irrespective of the value of k as the BH 
stage (2 = 00) would be arrived. Thus the Eddington-limited 
contracting phase actually becomes eternal. Since for pho- 
tons, K~i ~ 0.4 cm^/g, but, for neutrinos is smaller by an 
extremely large factor of ~ 10^*"^*, we will have Uv <^ u~f. 
Consequently, initial transition to the RRPSS phase may be 
dominated by huge t'-emission with a time scale Uv But as 
far as eventual secular RRPSS phase is concerned, it should 
be governed by photon time scale u~, because it is much eas- 
ier to maintain a L^d caused by photons than by neutrinos: 



1.3 



M 
IMn 



10^^(1-^2)-' erg/s 



(13) 



Somewhat similar thing happens for the formation of 
a hot neutron star from a proto neutron star: initial time 
scale of ~ 10s is dictated by huge j/-emission, while the hot 
NS cools for thousands of years by photon emission (Glen- 
denning 2000). For this era of quasi-stability by trapped 
photons, by using Eq.(13) into Eq.(12), it follows that the 
observed time scale of an RRPSS at a given 2 » 1 is given 
by 



1.5 X 10^^(1 -h 2) s w 4 X 10^(1 + 2) yr 



(14) 



For 2^1, the local energy density is almost en- 
tirely due to radiation and pairs (Mitra 2006a) so that 
p « aT^'/S where a is the radiation constant and T is the 
mean local temperature. Further since M = (4n/3) pR^ and 
R « 3{M/Mq) Km for 2 » 1 (see Eq.[l]), we obtain 



Q 4 \ 1/4 



600 



M 

M0 



-1/2 



MeV 



(15) 



Therefore, a RRPSS is an ultrarelativistic fireball of radi- 
ation and pairs interspersed with baryons much like the 
plasma in the very early universe (unless M is too high) . In 
particular, Eq.(15) shows that T 200 MeV for a lOM© 
RRPSS. Hence, the stellar mass RRPSSs could be in a 
Quark Gluon Plasma (QGP) phase. As of now, it is be- 
lieved that a bulk QGP phase existed only in the very early 
universe. But now we arrive, through a simple and straight 
forward analysis, at the exciting possibility that a bulk and 
ever lasting QGP phase may be existing within galaxies. 



4 ANALYTICAL & NUMERICAL SUPPORT 
FOR THIS SCENARIO 

The physical effect described hero cannot be obtained by 
any exact analytical solution of Einstein equations simply 
because the only exact solution of GR gravitational collapse 
is the Oppenheimer - Snyder one (Oppenheimer & Snyder 
1939): If one would consider the collapse of a homogeneous 
dust with p = 0, and yet assume finite initial density, p > 0, 
the fluid would appear to collapse to a singularity in a flash, 
T oc p~^^^, with no question of slowing down or bounce or 
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oscillation. However, physically, a strict p = EOS should 
correspond to a fluid mass of M = (Ivanov 2002, Mitra 
2009a) and thus, in a strict sense, a p = collapse should 
be eternal: r = oo. Note, when pressure gradient forces are 
included, even adiabatic GR collapse admits bounce and os- 
cillatory behavior where the fluid need not always plunge in- 
side its Schwarzschild radiu (Nariai 1967; Taub 1968; Bondi 
1969). Further Mansouri (1977) showed that a uniform den- 
sity sphere cannot undergo any adiabatic collapse at all if an 
equation of state would be assumed. Therefore, the effect 
described here here can be inferred by appropriate numeri- 
cal studies of radiative physical gravitational collapse or by 
generic physical studies, as carried out here. Indeed realistic 
gravitational collapse must involve not only pressure gra- 
dient but also dissipative processes and radiation emission 
(Mitra 2006a,b,c; Herrera & Santos 2004; Herrera, Frisco 
& Barreto 2006, Herrera, Frisco, & Ospino, 2006; Herrera, 
Frisco, Fuenmayor, & Troconis 2009). It is clear that dissi- 
pative processes might not only slow down but even stall the 
collapse by generating a quasistatic state (Herrera & San- 
tos 2004; Herrera, Frisco & Barreto 2006, Herrera, Frisco, 
& Ospino, 2006; Herrera, Frisco, Fuenmayor, & Troconis 
2009). And suppose, one is considering the collapse with a 
certain initial mass M = Mi and the RRPSS state is (first) 
formed at M = M* and z = z*. In order to find the val- 
ues of M* and z*, one must study the problem numerically 
by devising a scheme to incorporate the effect of gravita- 
tional radiation trapping and matter-radiation interaction. 
Obviously, the precise value of opacity k would be enor- 
mously different from what has been considered here. And 
same would be true for all other relevant physical param- 
eters. But at this juncture, we are not claiming to make 
any such detail study. Indeed, we arc interested only in a 
generic but physically valid picture. Implicitly, we are con- 
sidering here a modest range of Mi to relate it with asso- 
ciated rnimerical/analytical works (Herrera & Santos 2004; 
Herrera, Frisco & Barreto 2006, Herrera, Frisco, & Ospino, 
2006; Herrera, Frisco, Fuenmayor, & Troconis 2009). On the 
other hand, for extremely large values of M, there would be 
no z/-generation and one would be concerned with solely 7 
and e^,e~ processes. 



This generic/qualitative picture becomes strengthened 
by the fact Goswami & Joshi (2005) considered the possibil- 
ity that trapped surface formation may be avoided because 
of loss of mass by emission of radiation. Also, there is an 
exact solution of GR collapse which shows that the repul- 
sive effects of heat flow may prevent the formation of an 
event horizon (Banerjee, Chatterjee & Dadhich 2002). Fur- 
ther, Fayos & Torres (2008) have shown that in view of the 
emission of radiation, GR continued collapse may turn out 
to be singularity free where the entire mass-energy of the 
star may be radiated out. And recently, it has been found 
that, indeed, for radiative spherical collapse, trapped sur- 
faces are not formed at a finite value of AI (Mitra 2009b). 
Thus, there are considerable supports from both numerical 
and analytical studies of Einstein equations, about the basic 
feasibility of the picture presented here. 



5 SUMMARY & CLARIFICATIONS 

The generic, reason why continued radiative GR collapse 
should (first) result in a radiation supported quasistatic 
state or Eternally Collapsing Object(ECO), may be better 
appreciated by recalling Harrison: "When the contracting 
body reaches a radius 1.5 times the Schwarzschild radius, 
all rays emitted tangential to the surface are curved into 
circular orbits. This is the radius of the photon sphere. On 
further contraction, the emitted rays become more strongly 
deflected and many now fall back to the surface. Only the 
rays emitted within an exit cone can escape and this es- 
cape cone narrows as contraction continues. When the body 
reaches the Schwarzschild radius, the exit cone closes com- 
pletely and no light rays escape. Redshift and deflection con- 
spire to ensure that no radiation escapes from a black hole." 

Just like the case of an event horizon formation is inde- 
pendent of the density of the fluid, radiation trapping, the 
basic mechanism for the scenario presented here, too is in- 
dependent of density. It simply depends on the fact for con- 
tinued collapse, the fluid necessarily plunges within the pho- 
ton sphere and the effect of radiation trapping must increase 
dramatically with increasing, z. Essentially we have pursued 
this generic picture by recalling that there is another associ- 
ated generic concept - which is "Eddington Luminosity" . We 
showed that the rays falling back inside the contracting ob- 
ject should form an outward flux which would locally grow as 

+ whereas the definition of Eddington luminosity, in 
the increased gravity, would increase as ~ (I-I-2;). And since 
the former grows more rapidly by a factor of {l + z)/RM, it 
must catch up with the latter at some high z = z» ^ 1 as the 
journey to the BH stage involves march towards z 00 and 
A/ = (Mitra 2009a,b). True, there would be many non- 
generic effects like matter-radiation interaction, asphericity, 
rotation, magnetic field etc which would affect this collapse 
process. But all such effects would resist the free collapse sce- 
nario and hence they would reduce the value of a* . Thus this 
generic picture cannot be negated on the plea of complexi- 
ties of the GR collapse problem. In any case, publication of 
this paper may prompt researchers to incorporate this effect 
of gravitational trapping in numerical studies of radiative 
collapse. 

While the OS paper is the only truly exact study of con- 
tinued gravitational collapse, it is also the most unrealistic 
one because it assumes p = pr = q = even when the fluid 
would attain infinite density! And indeed Taub (1968) and 
Mansouri (1977) showed that a uniform density sphere can- 
not collapse at all if pressure would be incorporated! Con- 
sequently, in the past few years, several alternatives to BH 
Candidates have been proposed (Mazur & Mottola 2004; 
Chapline 2005). Unlike the OS case, such treatments are in- 
exact, and further they assume that the collapsing matter 
with positive pressure suddendly udergoes some unspecified 
phase transition to acquire negative pressure. Alternately, 
they assume that "Dark Energy" starts to play dominant 
role at local level by virtue of mysterios quantum processes. 
In contrast, though, our teatment too is inexact, we simply 
considered the natural GR effect that a contracting object 
must start trapping its own radiation once it is within its 
photon sphere. Note, one defines Event Horizon (EH) as a 
surface from which "nothing not even light can escape" . But 
the surface with Zc = — 1 is the precursor of an EH, 
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because outward movement of everything gets strongly in- 
hibited for z > Zc- Further, unless angular momentum is 
lost, everything may tend to move in closed circular orbits 
for z > Zc- However, there must be equal number of counter 
rotating orbits in order to conserve angular momentum. 

In such a case, all stresses must be tangential as radial 
stresses should almost vanish. In other words, the situation 
here may approach the idealized case of an "Einstein Clus- 
ter" (Einstein 1939). And it follows that, when stresses are 
completely tangential, there could indeed be static configu- 
rations with 2M/R 1 or 2 -i> oo Florides (1974). Thus, if 
one would wish, one might view a RRPSS/ECO as a (quasi) 
static configuration with z ^ 1. In such a case too, there 
is no upper limit on the value of M\ There is no denying 
the fact that we indeed observe compact objects, as mas- 
sive as 10^° Mq. It is also certain that compact objects of 
even 3 — 4M©s cannot be cold neutron stars. But it is fun- 
damentally impossible to prove that these objects are true 
BHs simply because, by definition, BH event horizons having 
z — oo cannot be directly detected (Abramowicz, Kluzniak 
& Lasota 2002). And as mentioned in the introduction, for 
sufficiently hot compact objects there is no upper mass limit 
even when stresses would be considered to be isotropic. As 
we found, for stellar mass cases, the quasistatic RRPSS fluid 
could be in a QGP state with T ~ 100 MeV. Interestingly, in 
the context of stellar mass BH formation, a recent study en- 
titled QCD Against Black Holes concluded that QCD phase 
transition would ensure that the collapsing object becomes 
a QGP fluid rather than a true BH (Royzen 2009). 

Note, there are many observational evidences (Robert- 
son & Leiter 2002; Robertson & Leiter, 2003; Robertson & 
Leiter 2004; Schild, Leiter, & Robertson 2006; Schild, Leiter 
& Robertson 2008) for the scenario described in this paper 
which show that the compact objects in some X-ray binaries 
or quasars could be RRPSSs with 2 S> 1 and strong intrin- 
sic magnetic moments rather than exact BHs with z — oo 
and no intrinsic magnetic moment. Also, by definition, it is 
impossible to claim that such compact objects are true BHs 
with exact EHs (Abramowicz, Kluzniak & Lasota 2002). 

Unlike Mazur & Mottola (2004) and Chapline (2005), 
the scenario considered here in no way denies that, mathe- 
matically, the final state of continued gravitational collapse 
is a BH. But, if one would ignore the phenomenon of trap- 
ping of own radiation due self-gravity as the body plunges 
into its photon sphere and also ignore the effect of pres- 
sure gradient or radiative transport, surely, one would find 
text book type prompt formation of BHs. However, the very 
fact that photon sphere is the precursor of an eventual EH, 
our scenario actually fills the missing gap between a photon 
sphere and a true EH during continued collapse. Accord- 
ingly, this paper may not be rejected on the assumption 
that it is in conflict with the basic mathematical notion of 
a BH. As the quasistatic hot RRPSSs would become more 
and more compact, they would asymptotically approach the 
ultimate state of spherical gravitational collapse; i.e., a BH 
with z = oo k M = (Mitra 2009a,b). 

Newtonian and Post Newtonian radiation pressure sup- 
ported stars, during their evolution, may be unstable to ra- 
dial oscillation (Chandrasekhar 1965). And there is a pre- 
liminary numerical computation which suggests that New- 
tonian supermassive stars may indeed collapse to form 
RRPSSs rather than true BHs (Cuesta, Salim & Santos 



2005). For Newtonian supermassive stars, even though. 
Pi- 2> Pg, p ^ pc^; in contrast the RRPSSs have p ~ 
(l/3)/9c^ Thus any study of Newtonian or Post Newtonian 
systems is not relevant for RRPSSs. Further, recently, it has 
been shown that, the Active Gravitational Mass Density of 
a quasistatic system is pg = p — ip/c? (Mitra 2010). Con- 
sequently, pg <^ p for RRPSSs and this prevents sudden 
rapid gravitational contraction. However, all RRPSSSs are 
assymptotically contracting towards the true BH state. 

To conclude, for continued collapse, we elaborated here 
on a most natural physical mechanism by which one may 
have GR version of Radiation Pressure Supported stars first 
conceived by Hoyle & Fowler way back in 1963. 
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